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Abstract 

We study the asymptotic behaviour of positive solutions of the Cauchy problem 
for the fast diffusion equation as t approaches the extinction time. We find a 
continuum of rates of convergence to a self-similar profile. These rates depend 
explicitly on the spatial decay rates of initial data. 
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1 Introduction 

We consider the Cauchy problem for the fast diffusion equation 

f u T = V • (u m - l Vu) , y £ W 1 , t G (0, T), 

\ u(y,0) = u (y) >o, j/er, 

where m < 1, n > 2 and T > 0. It is known that for m < m c := (n — 2)/n all 
solutions with initial data in some convenient space, like L p (M. n ) withp = n(l— m)/2, 
extinguish in finite time. We will always work in this range, m < m c , and consider 
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solutions which vanish in a finite time r = T. The purpose of this paper is to study 
the behaviour of such solutions near extinction. 

The situation is simpler and illustrative for our purposes in the exponent range 
m > m c that includes both fast (m < 1) and slow (m > 1) diffusion. Then solutions 
do not extinguish in finite time and the description of the asymptotic behaviour of 
the global-in-time solutions of (jl.ip as r — > oo for m > m c is a very active subject. 
For nonnegative data uq € L 1 (]R n ) it is proved (cf. |16| ) that all solutions converge 
up to scaling to one of the so-called Barenblatt solutions, precisely to the one with 
the same total mass, after proper rescaling of the solutions. For m c < m < 1 the 
Barenblatt solutions take the self-similar form 



Ud{v,t) : = (r + a) 



-fin 



D + 



13(1 - m) 



r + a 



1/3 



1 

1 — m 



(3:-- 



n \m c — m\ 



In other words, the Barenblatt solutions with exponent m > m c play the role of 
the Gaussian solution of the linear diffusion equation in describing the asymptotic 
behaviour of a very wide class of nonnegative solutions, i.e., those with integrable 
initial data, cf. [18]. The entropy method has allowed to provide rates for that 
convergence, cf. [6| \7\ l8j \W\ [TT] . 

This study has been extended recently to the behaviour near extinction for m < m c . 
Indeed, for m < m c and the problem posed in the whole space, the book [17] contains 
a general description of the phenomenon of extinction. It is explained there that the 
occurrence of extinction depends on the size of initial data, and also that different 
initial data may give rise to different extinction rates, even for the same extinction 
time; this may happen for all < m < m c . It is also proved in [T7J and references 
quoted there that the size of the initial data at infinity (the tail of uo) is very 
important in determining both the extinction time and the decay rates. 

Now, some natural questions arise, like: Are there solutions that continue the 
Barenblatt formulae for the extinction range m < m c ; if yes, how attractive are 
they? In fact, those solutions exist and have the self-similar form 



Ud,t(v,t) 



1 



R( T y 



D + 



13(1 -m) 



V 



R(r) 



1 

1 — m 



(1.2) 



where for m < m c we put R(t) := (T — r) ^, and 



n(l — m) 



n (m c — m) 



> 0. 



Many papers ( [21 El SI EQ], for example) are concerned with the stabilization as 
t —7- T of solutions of (jl.lj) . The general attractive character of the Barenblatt 
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solutions of the range m > m c is lost but still they have a basin of attraction 
formed by solutions with data that are close to the Barenblatt initial data in some 
norm. The study of such question is taken up in the papers [2j [3], which establish 
convergence with rates of the solutions of (jl.ip towards a unique attracting limit 
state in that family. More precisely, the decay rates of 

R(T) n (u(r,y)-U DiT (y,r)) 

as r — > T are discussed there (note that R(r) n U has size 1). 
The critical exponent 

n - 4 



has the property that the difference of two generalized Barenblatt solutions is inte- 
grable for m € (m*,m c ), while it is not integrable for m < m,. The exponent 
plays a very important role in the results of [21[3j|I]. The proofs of convergence with 
rates are based on the study of the decay in time of a certain relative entropy and a 
careful analysis of the linearized problem which leads to certain functional inequali- 
ties of Hardy-Poincare type. In particular, the basin of attraction of a function Ud,t 
in the range m < m* contains functions uq such that 

U Dl>T (;0) <u <U D2)T (-,0), \uo(x)-Ud,t(x,0)\ 6 L\R n ). 

We call this set the variational basin, and for this the entropy method gives precise 
decay rates (the variational rates). 

However we had the feeling that the basin of attraction is larger if we allow ourselves 
to get out of the variational framework, where the differences of solutions do not 
have finite relative entropy anymore. The analysis of this possibility was done in a 
very interesting limit case that occurs if we take D = in formula (jl.2j) , and we find 
the singular solution 

U , T {y, r) := h (T - TY' 2 \y\-^ K := (2(n - M ))"/ 2 , n. := . 

1 — m 

Its attractivity properties have been studied in [12] . confirming our guess in the form 
of a continuum of slow convergence rates for increasingly larger deviations at infinity 
with respect to the tail of Uo ; t, i-e., C|y| _At . 

The purpose of the present paper is to perform the construction of solutions in the 
basin of attraction of the generalized Barenblatt solutions Ud,t with both D, T > 
in the exponent range m < m*. Our solutions fall out of the already mentioned 
variational basin and consequently have a slower rate of convergence. 

To study the asymptotic profile as r — >■ T, it is convenient to rescale the flow and 
rewrite (jl.ip in self-similar variables by introducing the time-dependent change of 
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variables 

1, {R(t)\ 
t := — log — — and x : = 

with R as above, and the rescaled function 

v(x,t) :=R(T) n u(y,T). 



;i.3) 



;i.4) 



In these new variables, the generalized Barenblatt functions Ud,t{u^t) are trans- 
formed into generalized Barenblatt profiles Vd(x), which have the advantage of being 
stationary: 

V D {x):={D+\x\ 2 y 1 ^ 1 ~ m \ i£R". (1.5) 

The convergence theorems in the new variables take the form of stabilization to non- 
trivial equilibria. If u is a solution to (jl.ip , then v solves the rescaled fast diffusion 
equation 



v t = V 



,m— 1 



Vv) + juV • (xv), t>0, xe 



;i.6) 



which is a nonlinear Fokker-Planck equation. We put as initial condition Vq(x) := 
R(0)~ n uo(y), where x and y are related according to (|1 .3j) with r = 0. We have 
taken the precise form of this transformation from [2] . Note also that the factor fi in 
equation fjl .6|) can be eliminated by manipulating the change of variables, but then 
the expression of the Barenblatt solutions would contain new constants. 

By n > 2 and m < , we have (i + 2 < n so that the number 

n + u. + 2 



satisfies fi + 2 < I* < n. Note that as m — > m* we have \i 



n 



2 and L 



n. 
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m 

Figure 1. Values of Z* and fx + 2 as functions of m for n = 6 
Our main result is the following: 



Theorem 1.1. Let m < m* ; n > 2. Assume that c,D>0 and (i + 2 < I < I*. 
(i) U 

\v (x) - V D {x)\ < c\x\~ l , \x\>l, (1.7) 

and if 

< v {x) < Vg(x), x£R n 

for some 5 < D, then there exists C\ > such that the solution v of f 1 1 . 6 1) with the 
initial condition 

v(x,0) = v {x), x£W n , (1.8) 

satisfies 

sup \v{x,t) - V D {x)\ < Cie-^f 1 - 2 '^ 1 '', t > 0, (1.9) 
where (I — fj, — 2)(n — I) > 0. 
(ii) 

v (x) < V D {x) - c\x\~\ \x\ > 1, (1.10) 

and 

< v (x) < V D (x), x € M. n , 
then there exists C2 > such that the solution v of (jl.6p . (|1.8p satisfies 



sup (v D {x)-v(x,t)) >C 2 e- {l ->*- 2) ( n - l)t , t>0. (1.11) 

(iii) // 

u (x) > V D (x) + c\x\~ l , \x\>l, (1.12) 

and 

v (x)>V D {x), xeR n , 
then there exists C3 > such that the solution v of (|1.6|) . (|1.8p satisfies 



sup 



(v(z,t) - Vb(x)) > C7 3 e-f-"- 2) ( n - , J*, 4 > 0. (1.13) 



Remarks. (1) First of all, the result gives a sharp description of the basin of at- 
traction of generalized Barenblatt solutions for m < m*. For such m it was shown 
in [21 [3l [9] that the basin of attraction of a generalized Barenblatt solution Uu t 
contains all solutions corresponding to data which, besides being trapped between 
two Barenblatt profiles Ud ,t, Ud±,t for the same value of T, are integrable per- 
turbations of Ud,t- Theorem 11.11 (i) implies that if < uq(v) < Uspiy,®), 5 > 0, 
and 

\Mv) - u d ,t(v,o)\ < c\y\~\ \y\>i 
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for some c > and I > fj, + 2 then the corresponding solution is contained in the 

basin of attraction of Ud,t- Since m < m*, it follows that \i + 2 < n. Hence, non- 

integrable perturbations of £/dt may still yield convergence to Udt- The condition 

I > n + 2 is optimal since the difference of two Barenblatt solutions is of the order 
| y |-(M+2)_ 

(2) We have thus found a continuum of convergence rates which depend explicitly 
on the tail of initial data. The rate (I — fi — 2)(n — /) converges to zero as I — > \l + 2 
and to the maximum value 

as Z — > . Here a* is the rate found in [21 [3] for solutions emanating from integrable 
perturbations of Ud,t- This fastest rate is the best constant in a Hardy-Poincare 
inequality (see [3]). This best constant is also the bottom of the continuous spectrum 
of the linearization on a suitable weighted space (see [2j[3j[l0]). In the figure below, 
the gray area indicates the decay exponents that we obtain in Theorem II .11 for n = 6. 
The upper curve corresponds to the variational rate from [3]. 




-I -OJt -0,6 -0,4 -0.2 0,2 0.4 0.8 



Figure 2. 

(3) A continuum of rates of convergence to the singular Barenblatt profile £7o.r was 
found in |12j . For the existence of continua of convergence rates to positive steady 
states, zero and self-similar solutions of the Fujita equation 

u t = Au + u p , x 6 M n , 

see [13], [2] and [15], respectively. 

In Theorem II .1| the assertion (i) is no longer true if I > I*. In fact, we have the 
following result about the optimality of the range of I. 
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Theorem 1.2. Let m < m*. n > 2. Assume that D > ami 



< vq(x) < V D (x), 



x G 



or 



«o(a:) > V D (x), x G M n . 

Taen /or any e > 0, i/iere exists C £ > suc/i i/iai i/ie solution v of (jl.6p . (jl.8 
sate/ies 



sup 



> C e e- (a * +e)t , t > 0. 



(1.15) 



It follows from (|1.15j) that Theorem 2 (i) in [2] is optimal if m < m*, n > 2. The 
sharpness of the rate given by a* was discussed in [3] in terms of relative entropy 
which can be written as 



F[w] :-- 



1 



1 — m 



w 



1 



m 



-{w m - 1) 



Vffdx, 



w :- 



v 



The statement on the sharp rate in [3] says that a = a* is the best possible rate for 
which 

T[w(;t)]<T[w(;0)]e- at 

holds for all t > if Vb < f o < Vbj for some I?o > D > D\ > and vq — Vd 
is integrable. Theorem 11.21 implies that solutions starting from positive or negative 
perturbations of Vd cannot converge to Vd (in L°°) at exponential rates faster than 
e- a *K 

The paper is organised as follows. Section 2 contains some preliminaries, in Sec- 
tion 3 we derive an upper bound for the rate of convergence (Theorem ll.il (i)) and 
in Section 4 we show that the rate we found in Section 3 is optimal. We prove 
Theorem 11.21 in Section 5. Section 6 is devoted to a few comments. 



2 Preliminaries 



Let us consider the initial value problem 



/ TL — 1 

(r 2 + d) ( if rr H —¥r ) ~ /^r + Ctif = 0, V > 0, ^ ^ 



r 



\ p(0) = 1, <p r (0) = 0, 

with a parameter d > 0, which is related to the linearized eigenvalue problem at 
Vd(x). In order to study the behavior of the solution ip = ip d (r) of (|2.ip . we first 
consider properties of the differential operator 



CW := (r 2 + d) (w rr + ^— - W r ) - firW, 
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We define 

1(a) := n + , + 2-^{n-,-2Y-Aa > Q ^ 

for a € (0, a*]. We note that for a € (0, a*), 1(a) is the smaller root of the quadratic 
equation 

(I - pt - 2)(n - I) = a. (2.3) 

Lemma 2.1. Let a <E (0, a*] and W~(r) := (r 2 + d)~ fc/2 iwtift fc = 1(a) - fi - 2 > 0. 
T/ien i/ie inequality 

CW~ + aW - < 

/io/ds /or all r > 0. 
Proof. Since 

W~(r) = -Ax(r 2 + d) _fc/2_1 , 

W-(r) = jfe(ife + 2)r 2 (r 2 + d)"^ 2 " 2 - fc(r 2 + d)- fc/2_1 , 

we obtain 

= fc(]fe + 2)r 2 (r 2 + d)"^ 2 " 1 - £;(r 2 + d)~ fc / 2 - (n - l)k(r 2 + d)" k ' 2 

+ ll kr 2 (r 2 + d)~ k ' 2 - 1 
= (r 2 + d)-^ 2 - 1 (k(k + 2)r 2 - k(r 2 + d) - (n - l)k(r 2 + d) + juAx 2 ) 

= (r 2 + d)"^ 2 " 1 [(k 2 - (n - 2 - ^)A;)r 2 - nfcd) 

= (k 2 - (n - 2 - /i)fc)W~(r) - d(r 2 + d)- k/2 ~ 1 (k 2 - (n - 2 - fi)k - nk). 
Here, we have 

k 2 - (n - 2 - fi)k = (1(a) 2)(n - 1(a)) = -a 

by k = 1(a) — fi — 2, and 

k 2 - (n - 2 - fi)k - nk = l(a)(l(a) - fi - 2) > 
by 1(a) > [i + 2. Hence we obtain the desired inequality. □ 

Lemma 2.2. Let a G (0, a*) and W + (r) := r~ k — r~i with k = 1(a) — fi — 2 > 
and j = l((3) — fj, — 2. If /3 — a > is sufficiently small, then there exists ro > stic/z 
i/iaf f/ie inequality 

CW + + aiy+ > 

/ioWs /or r > ro. 



S 



Proof. By direct computation, we have 

Cr~ k = -ar~ k + 0(r^ k ~ 2 ), Cr~ j = -pr~ j + 0{r~ j - 2 ) as r -»■ oo. 
Since k<j<k + 2\{(3 — a > is small, we obtain 

£W + + alf + = (P — a)r~ j + o(r" J ) as r ->■ oo. 
This proves the lemma. □ 

The following lemma shall be frequently used throughout the sequel. 

Lemma 2.3. Let a G (0, a*). Then for every d > the solution ip = ip d (r) of (|2.ip 
is positive and decreasing in r G [0, oo). Moreover, the solution has the following 
properties: 

(i) There exist positive constants c and C such that 

cr -(K«)-M-2) < p«i( r ) < Cr .-(Ko ! )-/x-2) ) r > L (2.4) 

(ii) c/? d satisfies 

4 r (r) + ^-^(r) < 0, r > 0. 

(hi) satisfies 

0>4t4>-(H«)-^-2) , r>0. 

Proof. Define p(r) := r n ~ l (r 2 + d)~^ 1 ^ 2 . Then the equation in (|2.ip is rewritten as 
(p(r)^(r)) r + a(r 2 + d)-V(0/(^) = 0, r > 0. (2.5) 

Similarly, it follows from Lemma 12.11 that W~(r) := (r 2 + d)~ k l 2 with = /(a) — 
p — 2 > satisfies 

(p(r)W~(r)) P + a(r 2 + d)~ l p{r)W~{r) < 0, r > 0. (2.6) 

Multiplying (|2.5[) by VF~(r) and (|2 . 6[) by <p d (r), taking the difference and integrating 
it by parts on [0, r], we obtain 

p{r){ip d r {r)W-{r) - /(r)W"~(r)} > (2.7) 

as long as ip d (r) > 0. This implies that (p d (r)/W~ (r) is increasing in r, and hence 
(p d {r) > for all r > 0. Hence there exist Ci > such that C±W~(r) < ip d (r) for 
all r > 0. 
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Similarly, from Lemma 12.21 and (|2.5p , we have 



p(pfW + - <p d W+) < -a / (s 2 + dy 1 p(s)i P d (s)W + {s)ds. 

J r o Jr 

Here the integrand satisfies 

(s 2 + d)- 1 p(s)(p d (s)W + (s) > d(s 2 + d)- 1 p(s)W^(s)W + (s) 



(2.8) 



2 n-l-n -k -k 



(C 1 + o(l))s- 2 s 
(Ci + o(l))s n+ ^ +1 - 2 ' (a) 



s s 



as s — > oo. 



Since n + + 1 — 2Z(a) > — 1, the right-hand side of (|2.8p tends to — oo. This implies 
that (/? rf (r)/VF + (r) is decreasing for large r, and hence there exists C2 > such that 
(p d (r) < C2W + (r) for large r. This completes the proof of (i). 

Integrating (|2.5p on [0, r] and using the initial condition, we obtain 

p(r)<pt(r) = -a f\s 2 + d)' 1 p(s)y d {s)ds < 0, 
Jo 

which implies <p d (r) < for all r > 0, Hence by the equation in (|2.ip . (ii) is proved. 
Finally, by (|2.7|) . we have 

VrW ^ W r -(r) 



> 



-(1(a) -p -2) 



r 2 + d 



Thus (iii) is proved. 



□ 



Lemma 2.4. If a > a*, then the solution tp = (p d (r) of (12. ip changes sign. 

Proof. We will derive a contradiction by assuming that ip d (r) > for all r > 0. Set 
W*(r) := r~ fc with fe = l(a+) - fi — 2 = (n — p — 2)/2. Then we have 

£W* = -a*r~ fc + 0(r~ fc ~ 2 ) as r ->■ 00. 

This implies that for any fixed a € (0, a — a*), there exists r\ > such that 

(p(r)(W*(r)) r ) r + a(r 2 + d)' 1 p(r)W*(r) > a(r 2 + d) _1 p(r)W*( r ). r > r i- 

Then from (12.51). we obtain 



n 



p&fW* - /(W*) r )l < -a / (s 2 + d)~ Vs)/(s)T^(s)ds. (2.9) 
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Suppose here that {tp d {r) /W*{r)) r > holds for r > r±. Then there exists a 
constant c > such that tpd{ r ) > cr~ k for r > r\. In this case, the integrand in 
(1231) satisfies 

(s 2 + > (c + o(l)) S n ^" 2 s- 2fc = (c + o(l))s~ 1 as s -> oo. 

Hence the right-hand side of (|2.9p tends to — oo as r — > oo, which contradicts the 
assumption that (ip d (r)/W+(r)) r > for r ~>r\. 

Thus we may assume without loss of generality that (y d (^)/W / *(r)) r < at r = r±. 
Then by (|ZgD , we have 

p(r){^(r)W*(r) - /(r)(W*(r)) r } < -a / (s 2 + a!)-y S )/(s)^( S )ds (2.10) 

for r > n. Hence (</? rf (r)/W^(r)) r < for r > r\ so that 93 d (r)/W*(r) converges 
to a nonnegative constant. This implies that there exists a sequence {r^} such that 
r, — > oo as i — > oo and (<^ rf (r)/W^(r)) r — > along the sequence. Then the left-hand 
side of (|2.10p tends to along the sequence, since p(r) and W*{r)~ 2 are of the same 
order as r — >■ oo, whereas the right-hand of is negative and decreasing in r > r±. 
This contradiction shows that <p d (r) must change sign at some finite r. □ 

Let V be the differential operator 

(w m ~ l w r ) r H w m ~ l w r j - Lirw r - fmw, V > 0, t > 0. (2.11) 

In the following lemma we evaluate Vw for a suitable auxiliary function w which 
will be used frequently in the sequel. 

Lemma 2.5. Let D > 0, and y : [0, oo) — > K. and ip : (0, oo) — > (0, oo) be any 
smooth functions. Then the function 

w(r,t) := (r 2 + D + y(t)ip(r)j , r > 0, t > 0, 

sato/ies 

Pi0 = fyW(V 2 + £> + y(t)^(r)) " (M+2)/2 ^ D [y(t)]^, r>0, t > 0, (2.12) 

A fl [y(t)]V := (r 2 + £>) (VVr + — -A) ~ ~ ip 

-y(t)| + + ^(W) 2 } (2.13) 

/or r > and i > 0. 
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Proof. We compute 

w t = --[r + D + y(t) tp(r)J y'tp 

and 

w r = "f( r + D + y(t)^(r)j (2r + yW) 

and using the relation m = (// — 2)//i we find that 

w m-l Wr = _Ji_ ( (r 2 + D + ^)-(M-2)/2 



r 



= -|(r 2 + £> + y(i)^(r)) ^r + yW) 
as well as 

( W m - V)r = ^{r 2 + D + y(t)^(r))" (At+2)/2 (2r + y?/v) 2 

- |(r 2 + Z) + y(i)V>(0) ~ W2 (2 + # rr ) 
for r > and f > 0. Therefore, 
Vw = |(r 2 + D + y{t)i>{rj) (m+2)/2 | _ - |(2r + yW) 2 

+ (r 2 + D + yV')(2 + yVrr) 

+ (r 2 + D + yV)(2r + yVv) 

r 

+ /xr(2r + yVv) — 2n(r 2 + D + yVO} 

= |(V 2 + D + yWV'CO)"^ 72 ! - yV - 2^r 2 - 2//yr?/v - |y 2 (Vv) 2 

+ 2(r 2 + D + yip) + y(r 2 + d + yVOVVr 
+ 2(n- l)(r 2 +L> + yV>) 

+ y(r 2 + D + yV)- — -^r + 2 A«" 2 + 
r 

-2n(r 2 + J D + y^)}, 
which after obvious simplifications yields (I2.12p . □ 

3 Upper bound for the convergence rate 

The goal of this section is to provide an upper estimate for the distance between the 
solution v of (jl,6p and Vd under the assumption that |t>o(V) — Vo( r )| < cr~ l for all 
r > 1 and some c > and / € (fi + 2, with given in Theorem II. li 
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3.1 Estimating v(-,t) — Vd from above 



We first consider the positive part of the distance, that is, we seek an upper bound 
for sup r > (f (r, t) — Vb(r)) for large t. This will be accomplished along a three- 
step procedure: In a first step, in Lemma 13.21 we shall establish that this term will 
be sufficiently small, but at a rate possibly smaller than the desired one. Next, 
in Lemma 13.31 we shall make sure that the spatial decay rate of the 'positive tail' 
(v(-,t) — Vd)+ is essentially preserved during the evolution. Finally, the proof of 
Lemma 13.41 will be based upon the observation that if at a certain time the solution 
is close enough to Vd and its positive tail decays sufficiently fast in space, then 
(v(-,t) — Vd)+ will decay in time at the desired rate. 

In order to avoid redundant arguments, we start with the following lemma that will 
be used twice in Lemma 13.21 and in Lemma [ 



Lemma 3.1. Let D > 0,5 G (0,D),B > and a G (0, a*) with a+ given by (|1.14|) . 
and suppose that ip s solves (|2.ip with d = 5. Moreover, assume that rj is a real 
number satisfying 

r } <a-^-^(l(a)-^-2]\ (3.1) 



25 

where 1(a) is defined by (|2.2p . Then the function v defined by 

v(r,t) :=mm{v 5 (r) , (r 2 + D - 5e _r? V (0) ^ J > r > 0, t > 0, (3.2) 
satisfies Vv > in the Nagumo sense in (0, oo) x (0, oo). 



Proof. Since is a solution of (jl.6p , we only need to make sure that Vv > is valid 
at each point from the set 



S 



:= |(r,t) G (0,oo) x (0, oo) | v{r,t) < ^(r)}. 



In view of Lemma [231 this amounts to showing that at such points Ad[ — Be ^j^ 5 < 
holds, where 

Aol-Be-*]^ = (r 2 + D)(4 r + ^rf)-»rtp s r + r,<p s 

-vtf _,J(,J ■ n ~ 1 .J\ ■ Pt.JSV 



+Be*\ - + — ) + Z(tfyy (3-3) 

To see this, we first note that the nonlinear terms —(p s ((pf r + ^j^tPr) an d ((fir) 2 
in (|3.3p are both positive on (0, oo) by Lemma 12.31 In order to estimate them 
appropriately, we make use of the fact that whenever v(r,t) < Vg(r), due to (j3.2[) 
we know that 

r 2 + 5 = V s {r)- 2 ^ < v(r, t)" 2 ^ = r 2 + D - Be^/(r) 
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and hence 

D-S 



Be~ vt < 



(p 5 (r) ' 
Using this in (|3,3p yields 

A D [-Be-^ 5 < (r 2 + D)(tp s „ + ^¥>?) - »rrf + T)<p s 
-(D-d)[ ip rr H ip r + 



(r 2 + 5) <p% + 4 - urtp* + w 5 + H 



r ttj r- rr ■ -,r ■ 2(fiS 

at such points. According to the ODE in (|2.ip . we thus obtain that 

2</r 



2<p s 



is satisfied at all points in 5. Now from Lemma 12.31 (hi) we know that 

<5\ 2 „ o 

^ r ' < f /(a) - /i - 2 



ip 5 J ~ \ K ' h J (r 2 + 5) 2 

/ \ 2 1 

< ( Z(a ) _ ^ _ 2 j 



r 2 + <5 



1 / \ 2 
< - [1(a) -fi-2) for all r > 0, 

which inserted into (|3.4p yields 

>M-Be-"V <{-(«-»/) + M(D 2 ~ 5) (/(a) - /i - 2) 2 |/ in 5. 

This immediately proves that under the assumption (|3.ip . *4_d[ — -Be -77 *]^" 5 < holds 
throughout S. □ 



As announced, our three-step procedure is launched by deriving an estimate from 
above. 

Lemma 3.2. Suppose that there exist 5 > 0, D > 8, c > and I £ (/i + 2, £*) suc/i 

< uo(r) < F 5 (r) /or all r > (3.5) 
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and 

vo( r ) < Vd(t) + cr~ l for all r > 1. (3-6) 

T/ien tZiere exist C > 0, ry > and Iq G (/j + 2, Z] sitc/i i/iai £/ie solution v of U.6\) 
satisfies 

v(r,t) < V D {r) + Ce~ vt (r + for all r > and t > 0. (3.7) 

Proof Since our assumption m < (n — 4)/(n — 2) implies /i + 2 < n, we can pick 
some Iq > fx + 2 close enough to /i + 2 such that lo < I, Iq < I* and 

n — to — o) 

Then a := (Iq — [i — 2){n — Iq) is positive, so that we can pick some rj € (0, a/2), 
and by (|2.2p Zo is given by 

Z = Z(a). (3.9) 
Now from the hypothesis (|3.6p and the convexity of < z i-> z~ 2 /^ we obtain 



wo 



( r )-2/A* > (yb(r)+cr _i ) 



-2/n 



= r 2 + D _ ^) ( r 2 + D) (M+2)/2 r -Z for an r > 1; 

which combined with the positivity of Vd and fj2.4j) easily yields that 

Vp 2/f *(r) >r 2 + D - Bip S (r) for all r > (3.10) 

holds with some suitably large B > 0, where ip s denotes the solution of (|2.ip with 
d = 5. In view of (|3.5p . this implies that the function u defined by (13. 2p dominates 
w at i = 0. Since moreover 

'Z(a)-/x-2 = Z - ^ - 2 = — — — < 



-i 



26 V v ' *~ ) 28 V u ^ ) 2(n- 1 )5 ~ 2 

by (|3.9p and (|3.8p , Lemma 13.11 yields that v is a subsolution of (jl.6p and hence 
we have i> < u by comparison. Now (|3.7p is a straightforward consequence of the 
definition of v: We let ci > be such that (1 - z)~^ /2 <l + c x z for all z G [0, £/£>] 
and observe that since < we have 

Be-^(p 8 (r) B 
r 2 + D ~ D 
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for all r > and t > 0. Hence, 

v(r,t) < v(r,t) 



< [r 2 + D-Be- Tlt v 5 {r)^ 



(r 2 + D) ^/2 L 



r 



-M/2 



r 2 + D 

< {r 2 + D)~» 12 + c l Be- r > t {r 2 + D)-^ +2)l2 v\r) (3.11) 



for all r > and i > 0. By means of the second inequality in (|2.4p we easily arrive 
at (J373). □ 

We next assert that the decay of (v (•, t) — Vd)+ does not change in time significantly. 

Lemma 3.3. Under the assumption of Lemma \3.S\ there exist C > and k > 
suc/i the solution v of (|1.6|) satisfies 

v(r, t) < V D {r) + Ce K *(r + 1)"' /or all r > and t > 0. (3.12) 

Proof. Since I £ (/" + 2, £*), the number a := (/ — fi — 2){n — I) is positive and satisfies 
a < a* with a* as in (jl,14j) . and moreover the number 1(a) in (|2.2p coincides with 
L We now choose positive constants n and such that 

u(D-5){l-y-2) 2 

k> — a (3.13) 

and 

v~ 2/fl (r) >r 2 + D - Btp 5 (r) for all r > 0, (3.14) 

where the latter is possible because of (13.6|) and (12. 4h (cf. the argument leading to 
fl3~T0]) ). Thus, defining t> by (J3T2]) with 77 := -k, (pTT4"P states that u (r) < v(r,0) 
for all r > 0, whereas (|3.13p ensures that (13. If) holds. Therefore, Lemma 13.11 in 
conjunction with the comparison principle guarantees that v < v in [0, oo) x [0, oo), 
whence (|3.12p easily results from the definition of v (see (|3.10p ) . □ 



We are now in the position to prove the main result of this section. 
Lemma 3.4. Under the assumption of Lemma \3. t A the solution v of (|1.6|) satisfies 

v(r, t) < V D (r) + Ce- {l -^- 2)(n - l)t (r + l)"^" 2 ) f or a ll r > and t > (3.15) 

with some C > 0. In particular, v satisfies 

sup (v(r, t) - V D (r)) < Ce^ 1 "^ 2 ^-^ for all t > 0. (3.16) 
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Proof. Again, at := (I — fi — 2)(n — I) belongs to (0,a*) and is such that 1(a) = I. 
Now taking any small 5 > such that 

X ^ (n + n-l)D 

0< n + n-l + 1(l-n-2y [6An 

from Lemma 13.21 we infer the existence of some large to > such that 

v(r,to) < V D -s{r) for all r > 0, (3.18) 

because the difference Vz)_s(r) — Voir) decays with the order r~^ +2 \ Then in view 
of Lemma 13.31 it is possible to fix a large number c\ > such that 

v(r, t ) < V D (r) + cir~ l for all r > 1. (3.19) 

We finally choose B appropriately large satisfying 

B>e ato (D + l) (3.20) 



and 



and let 



B > ^l e «*o (jD + 1) (^+2)/2 (3 _ 21) 



v(r,t) := min|vb-«(r) , (r 2 +I}-£e _a V~ ( ^ _2) ) + M/2 J, r > 0, t > 0. (3.22) 

We first claim that v > v holds for all r > at t = to- Indeed, if r < 1 then 

r 2 + D - Be" Qi °r- (/ ~^ 2) < 1 + D - Be~ at ° < 

by (|3T20|) . hence t>(r,t ) = Vo-^r) > u(r,*o) by (ETTKD . On the other hand, if r > 1 
then by convexity of < z i— y z~ 2 l^ and (|3.19p . 



^(r,t )~ 2/M > (v D {r) + cr- 1 ^ 
= V D {r)- 2 l»(l + 



Vb(r) 

> Vb(r)- 2 /" - —V D (r)- 2 l»- x r- 1 

9 2c fr 2 +D\(m+2)/2 - . 
= r 2 +L> ( 5— ] r -('-M-2) 

>r 2 + D — — (D + i)0*+2)/2 r -d-M-2) 
A* 

>r 2 + D- Be~ ato r~^^~ 2) for all r > 1 
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in view of (|3.2ip . 

Having thereby asserted the desired ordering at t = to, let us next make sure that 
v is a supersolution of the PDE in (II. 6j) for r > and t > to- To achieve this, since 
Vds solves (jl.6p . we only need to check that Ad[— Be~ at ]ip < in 5 is valid with 

ijj(r) ■= r-C-^" 2 ) and S := {(r,t) G (0, oo) x (t ,oo) | u(r,t) < Vd-<5(V)}; where 

+ fle-^j - ^ rr + ^V) + 
Differentiating tp, we immediately obtain 

A D [-Be~ at ]iP = -(r 2 + D)(l - fi - 2)(n + fi - 0r _M 
+ fi(l -fi- 2)r- (l -^ + ar"^^- 2 ) 
+ Be~ at \r-^ l -^- 2 \l -fi- 2)(n + fi - Z)r"^ 

+ |((/-^-2) 2 r"^)) 2 } 

= -(Z - // - 2)(n + [i - l){r 2 + D)r^ l ~^ l) 
+ (^(/- / u-2) + a)r"^- 2 ) 

+ 5e- at (/ - fi - 2)(n + /i - I + |(Z - /i - 2))r~ 2 ^- 1 ), 
so that by definition of a we see that 

A D [-Be~ at ]tp = -(l-n-2)(n + n- Z)I> _M 

at (l — n — 0\( n -i- ll — 7-1- —(1 — il — 0\\ r -2(l-iM-X) 



+ Be' a \l - fi - 2)(n + fi - I + ^{l - fi - 2))r~ 

for r > and t > 0. Now if (r, t) G S then, by (pT22l . 

5e~ Q< < 5r'-^" 2 , 
and therefore at these points we have 

A D [-Be- at }4> <-{l-fi- 2)(n + fi - l)Dr~^ 



+ 5(l-fi-2)(n + fi-l + ^{l-fi-2^j r-Q-ti 
(Z-/i-2){- (n + n-l)D 
+ <f(n + /x-Z + |(Z-/x-2))}i 



r -(J-m) 



< 



because of our smallness condition A3. lTf) on <5. Consequently, by comparison we 
conclude that v < v for r > and i > £o 5 an d hence (|3. 15 j) and (|3.16p easily result 
from the definition (|3.22p of v. □ 
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3.2 Estimating v(-,t) — Vd from below 

A corresponding statement on convergence from below can be obtained much more 
directly by a one-step comparison argument. 

Lemma 3.5. Suppose that vq > on [0, oo), and that there exist D > 0, c > and 
I G (/i + 2, 1+) such that 

i>o(r) > V D {r) - cr~ l for all r > 0. (3.23) 

Then there exists C > such that the solution v of (|1.6|) satisfies 

v(r,t) >V D {r)-Ce- {l -v-^ n - l)t {r + l)- 1 for all r > and t > (3.24) 

and £/ms 

sup (Vb(r) - v(r, t)) < Ce-V-t*- 2 ^- 1 ^ for all t > 0. 

Proof. Since I > fj, + 2 > fj,, it is possible to fix ro > 1 such that with c as in (|3.23D 
we have 

cr -l < -V D {r) = -(r 2 + D)-»/ 2 for all r > r . (3.25) 

Then according to Lemma 12.31 there exists c\ > such that the solution of (|2.1|) 
with d = D satisfies 

ip D (r) > cir~ (i ~ At " 2) for all r > r , (3.26) 
and using that vq > on [0, oo) we can pick C2 > such that 

«o(r) > c 2 for all r € [0, r ]. (3.27) 
We next let C3 > be such that 

(1 - z)~l < 1 + c 3 z for all z G [0, 1/2] (3.28) 
and finally fix B > large enough satisfying 

B> (^V(ro))" 1 (3.29) 

as well as 

7(g + 1) W 

ij > , 7 := CC3. (3.30) 

ci 

Then, again writing a := (I — fi — 2)(n — Z) > 0, from Lemma 12.51 we know that 
t;(r, t) := (r 2 + D + Be~ at ip D (r)J ^ , r > 0, t > 0, 
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defines a subsolution of (jl.6p if and only if An\Be ]ip < 0, where 
A D [Be- at ]^ D = (r 2 + D)L° + ^— ^) - /^f + c^ D 



Be"^ - + ^rV) + \ \, r>0,t>0. 

Recalling that Lemma [231 yields (p^ r + ^r-tp® < on (0, oo), in view of (|2.12p we 
directly see that indeed 

A D [Be- at ]<p D < (r 2 + D) (<p° + ^— V) - ^ + a<fP = 



for all r > and i > 0. In order to use v as a comparison function, we next claim 
that 

v(r, 0) < v Q (r) for all r > 0. (3.31) 
In fact, for small r, (f3T27|) and (pl29|) tell us that 

, M )^ = l fr2 + I) + BA)r ,/ 2 



u (r) c 2 c 2 

1 / n X-A 4 / 2 

<-£<Ar <1, rG[0,r ], 

c 2 V / 

because </? D decreases. For large r, we may use (|3.25p and (I3.28P to estimate 

voir)' 2 /" < (Vb(r) - cr~ i y 2/fl 

< V D (r)- 2 ^ + 7 V D (r)- 2 ^- 1 r- i 

2 



2 /H + D\(M+2)/2 

r 2 + L» + 7( ^ — ) r 



< r 2 + D + 7 (D + l)0*+2)/2 r -(I-M-2) for all r > rQj 

whereas 

tj(r, O)' 2 ^ =r 2 + D + B(p D (r) 

>r 2 + D + ci-Br _ ^ _M_2 ) for all r > r 

by (I3.26P . Therefore (I3.30|) completes the proof of (|3.31|) , which in turn implies that 
v > v holds for all r > and t > by comparison. As a consequence we infer that 
if we pick C4 > small such that (1 + z)~^/ 2 < 1 — c^z for all z E [0, B/D] then 

v(r, t) < (r 2 + D)-»/ 2 - c A Be~ at {r 2 + D)~^+ 2 )/ V \r), r > 0, t>0, 

which immediately implies (I3.24p for some C > upon invoking (|2.4p . □ 
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4 Optimality 



We proceed to show the optimality of the convergence rates that we have found in 
the previous section. 

Let us first make sure that convergence from below does not occur at a rate faster 
than asserted by Lemma 13.51 

Lemma 4.1. Suppose that there exist D > 0, c > and I G (fi + 2, /*) such that 

v (r) < V D (r) for all r > (4.1) 

and 

Vo(r) < Vo{r) - cr~ l for all r > 1. (4.2) 
Then there exists C > such that the solution v of (jl.6p satisfies 

v(r, t) < V D (r) - Ce~ {l ~ ll ~' 2){n ~ l)t {r + for all r > and t > 0. (4.3) 

In particular, v satisfies 

sup (v D (r) - v(r, t)) > c e -('-M-2)(n-«)t f or aU t > . (4.4) 

Proof. In view of (|4.ip and ()4.2p it is possible to fix a number E > D close enough 
to L> such that 

and 

u (r) < Vk(r) for all r < 1, (4.6) 

where again a := (I — \i — 2)(n — I) > 0. Applying Lemma [2. 3 \ we obtain a constant 
c\ > such that the solution of (|2.1|) with parameter d = E + 1 satisfies 

^ +1 (r) < Cl r-(^- 2 ) for all r > 1, (4.7) 

and finally we pick B > so small that 

Then 

^r,t) :=max!^V E {r), (r 2 + D + Be- at <p E+1 (rfj ^j, r > 0, t > 0, (4.9) 
satisfies 

w(r, 0) > Ve(^) > ^o( r ) fo r all r < 1 
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by (USD, whereas (|3~2j) . (J1771) . (|3~8j) and the convexity of < z \-+ z 2 ^ ensure that 

U(r,0)" 2//1 -v (r)- 2/ » < (r 2 + D + By E+l (r)} - (v D {r) - cr - ')" 27 " 

< (V 2 +,D + S^ £;+1 (r; 

r 2 + £)+ 2 f( r 2 +jD) -(M+2)/2 r -i 



r; „ s>i 2c/ r 2 \ 0+2 )/2 , 
fi \r 2 + DJ 
< for all r > 1, 



because r 2 /(r 2 + -D) < 1/(D + 1) for such r. We thus conclude that v > v at t = 0, 
and next claim that > for r > and i > 0. According to Lemma 12.51 and the 
fact that Ve is a solution of (jl.6p . we need to check that for ip := tp E+1 , 

A D [Be- at ]if; = (r 2 + D) (^ rr + ^r~Vv) - A""Vv + «V> 



Be~ at { -1p(lp r r + — -A) + ^A 



A + 7T-'- 

r J 2 

is nonnegative whenever v(r, t) > Ve(t), that is, when Be~ at ip(r) < E — D. At such 
points, using (|2.ip we have 

-oti./. ^ a_2 ■ — t • n-\ 



+ (E-D)[ lfj rr + —A 

n — 1 



A D [Be at ]ip > {r 2 + D) [ Ar + —A) ~ ^A + on/; 

n-l L \ n{E — D)A^ 

At h — - — A) - v r A + Otlp 

il{E - D)Ar 



2ip 



fxT a fi(E - D)A r 



A + 9 , „ , > 



r 2 +E + V r 2 + E+V 2i/j 



a fi(E - D)ip. 



r 2 + E+l T 2ip 
since ip decreases. Thus, in view of Lemma 12.31 and (|4.5|) we infer that 

H(l - H-2) 2 (E - D) r 2 \ V 



A D [Be- at \^ > < a 



> <a 

> 



r 2 + E+l\r 2 + E+l 
H(l -fi- 2f{E -D)\ ip 

r 2 +E + l 
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is valid whenever v > Ve, which shows that indeed v is a supersolution of (jl.6p . so 
that the comparison principle states that v < v in [0, oo) x [0, oo). Since 

r 2 + D + Be~ at ip{r) < r 2 + D + Be~ at < r 2 + E for all r > 

when t is sufficiently large, using (|4.9p we easily derive (|4.3p and (|4.4j> . □ 

A similar result on the rate of convergence from above can be obtained in a slightly 
easier manner. 

Lemma 4.2. Assume that 

i>o(r) > V D {r) for all r > (4.10) 

and 

Vo(r) > V D (r) + cr~ l for all r > 1 (4.11) 
are valid with some D > 0, c > and Z € (/i + 2, I*). Then there exists C > swc/i 

u(r, t) > V D (r) + (7e - ( ,- ' 4 - 2 ^ n - ,)t (r + for all r > and i > 0, (4.12) 

and in particular 

sup (v(r, t) - V D (r)) > c e -(^-2)(n-fj< j or a // f > . (4.13) 

Proof. According to (I4,10p and (14. lip , it is easy to see that for some c\ > we have 
Vo{r)- 2/ » < V D {r)~ 2 ^ - Cl (r 2 + £))-(l-/*-2)/2 for all r > _ 

In view of Lemma 12. 3| this means that if B > is small enough then taking tp D 
from (]2.ip with d = D, we have 

^o(0~ 2/m < Vd(?")~ 2/m - Btp D (r) for all r > 0. 
Hence, the function v defined by 

v(r, t) := (r 2 + D - Be~ at y D (r)\ ^ , r > 0, i > 0, 
again with a := (I — fi — 2) (n — /) > 0, does not exceed u at i = 0. Since 

A fl [-5e- Qt ]^ = (r 2 + D) (y£ + ^V D ) " M^f + <*<p D 

+ 5 C -{ - + ^V) + | (^) 2 } 
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satisfies 

A D [-Be- at ]<p D > (r 2 + D)^ r + - l*r<p? + ay D = 

for all r > and t > by Lemma 12. 3[ we conclude using Lemma 12.51 and the 
comparison principle that v > v for r > and t > 0, which easily results in (|4.12p 
and (14131) . □ 



5 Universal lower bound for the convergence rate 

Proof of Theorem \1.2l Let us first assume that < vq < Vr> in W 1 for some D > 0. 
Then given e > 0, we write a := a± + e and let 99 = ip d (r) denote the solution of (|2.1|) 
corresponding to this value of a and d := D + 1. Then since a > a+, Lemma 12.41 
says that there exists ro > such that (p(r) > for r G [0,ro) and y?(ro) = 0. 
Moreover, from (|2.ip it is immediately clear that y? r < in (0, ro), because at each 
point in (0, ro) where ip r vanishes we must have (p rr < by positivity of a and if. 
Since (p r {ro) cannot vanish in view of an ODE uniqueness argument, it follows that 
<p r < even in (0,ro]. This in turn entails that 

n — 1 finp r — oup 

ip rr H ip r = — — - — < -ci in (0, r ) (5.1) 

r r z + a 

is valid for some c\ > 0, so that by smoothness of ip, 

C2 ■= { Cl r (5.2) 

SUp re(0 ,r ) yP\Vrr + + %<Pr\ 

is positive. 

Now since vq(x) < Vd(x) = (|x| + D)~z in M n and vq is continuous, it is possible 
to fix i/Q G (0, C2] small enough such that 



vq(x) < l\x\ +D + y f{\x\)j for |sc| < r , (5.3) 

whence the function 

v(r,t):= (r 2 + D + y{t)^{r)} \ (r,t) G [0,r ] X [0, 00), 

with 

y(t) := y e~ Qi , i > 0, 
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satisfies 

Vq(x) <v(\x\,0) for \x\ < r . (5-4) 

Since clearly vq < Vb i n implies v(x,t) < Vd(x) in M. n x (0, oo) by comparison, 
we also have 

v(x,t) < Vd(x) = T>(|x|,t) whenever \x\ = r and i > 0. (5-5) 
Furthermore, recalling the definition (|2.1ip of V, by Lemma 12.51 we have 

^= |y(*)(^ 2 +-D + y(*Mr)) 2 A D [y(t)]ip for r € (0,r ) and t > 0, (5.6) 
where 

ip rr H — VrJ - A""¥V - — V 

™ - 1 \ , A 1 2~l 



/ n — 1 \ /x 2 
¥\Wr H — J + 9^ 



(D — d) ( (/? rr + - y r ^ + ( - - — ajtp 
n-1 \ u 2 ' 



according to (|2. 1 j) . Since ci = Z? + 1 and y'/y = —a, we thus have 

•AD[y(t)]<P = - (Wrr + ~~~^) ~ y( t ){ ~ <p(<Prr + "~~"^) + f^r} 



f n — 1 ^21 

> Ci - yo|^ Vrr H — +2"^ J 

> for r e (0, r ) and t > 



because of (|5.ip . f)5.2j) and the fact that j/o < C2. In view of (|5.4p - (|5.6p and a 
comparison argument this implies that 

v(x,t) <v(\x\,t) for \x\ < ro and i > 0, 
which in particular means that 

Vb(0) - u(0, t) > D~2 - (D + y e~ at )~2 for all t > 0. 
Now an application of the elementary inequality 

{1 + z)~ 2 <\-^-(\ + zqY^z, 0<2;<zo<oo, 
to z := ZQe~ at and zo := 2/o/-^ shows that 

M+2 



Vb(0) - u(0, t) > ir£|(l + g) 2 e - Q * for all* >0 
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and thereby establishes (|1.15p in the case when < vq < Vr>. 
The proof in the case vq > Vd can be carried out along the same lines. Indeed, 

i£ 

v(r,t):= (r 2 + D-y(t)p(rf} \ (r, t) £ [0, r ] x [0, oo), 

becomes a subsolution of (jl.6p satisfying v(x,t) > v(\x\,t) in R n x (0, oo), provided 
that ro and <p are chosen as above and y(t) := yoe~ at , where a is as before and 
jjo > is small enough fulfilling 

Vo(x) > (\x\ 2 + D — yQip^xD^j 2 for |x| < ro 
and yo < min{c2, D} with C2 as given by (|5.2p . In conclusion, in this case we obtain 

v(0,t)-V D (0) > (D-y e- at Y* - D~$ 

> D'^-e- at forallt>0 
2 

by convexity of [0, 1) 9 z i— > (1 + z)~% . This again proves (|1.15p . □ 

6 Comments 

We have not discussed the problem of constructing solutions in the non-variational 
basin for m > m*. For m > m*, the existence of such a basin seems unlikely but 
for m = m* we expect that the variational basin studied in [3] represents just a part 
of the whole basin of attraction and the variational rate found in [3] is the fastest 
possible rate of convergence. The construction of solutions with slower rates for 
m = m* seems challenging and it will become a theme of future investigations. 

The questions of extinction behaviour are different when the problem is posed in 
a bounded domain. Then the extinction rate for bounded solutions is universal, of 
the form ||u(-,t)||oo = 0((T - r) 1 ^ 1 "" 1 )) when m > m s := (n-2)/(n + 2), cf. p], 
but the question is more complicated when m < m s . Convergence rates for m near 
1 have been recently obtained in [5]. 
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